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Abstract 

Stochastic gradient descent (SGD) is a ubiquitous algorithm for a variety of machine learning problems. Re¬ 
searchers and industry have developed several techniques to optimize SGD’s runtime performance, including asyn¬ 
chronous execution and reduced precision. Our main result is a martingale-based analysis that enables us to capture 
the rich noise models that may arise from such techniques. Specifically, we use our new analysis in three ways: (1) we 
derive convergence rates for the convex case (HOGWILD!) with relaxed assumptions on the sparsity of the problem; 
(2) we analyze asynchronous SGD algorithms for non-convex matrix problems including matrix completion; and (3) 
we design and analyze an asynchronous SGD algorithm, called Buckwild!, that uses lower-precision arithmetic. 
We show experimentally that our algorithms run efficiently for a variety of problems on modern hardware. 


1 Introduction 


Many problems in machine learning can be written as a stochastic optimization problem 

minimize E[/(x)] over x £ IR n , 

where / is a random objective function. One popular method to solve this is with stochastic gradient descent (SGD), 
an iterative method which, at each timestep t, chooses a random objective sample f t and updates 


x t+ i =x t - aVf t (x t ), 


(1) 


where a is the step size. For most problems, this update step is easy to compute, and perhaps because of this SGD 
is a ubiquitous algorithm with a wide range of applications in machine learning [7|], including neural network back- 
propagation 1128.13, 11311 . recommendation systems |8|,|l9|], and optimization |2Cf]. For non-convex problems, SGD is 
popular—in particular, it is widely used in deep learning—but its success is poorly understood theoretically. 

Given SGD’s success in industry, practitioners have developed methods to speed up its computation. One popular 
method to speed up SGD and related algorithms is using asynchronous execution. In an asynchronous algorithm, 
such as HOGWfLD! [17], multiple threads run an update rule such as Equation Q] in parallel without locks. HOG- 
WfLD! and other lock-free algorithms have been applied to a variety of uses, including PageRank approximations 
(FrogWild! [JLSD, deep learning (Dogwild ! H) and recommender systems [[24]. Many asynchronous versions of 
other stochastic algorithms have been individually analyzed, such as stochastic coordinate descent (SGD) 01 
and accelerated parallel proximal coordinate descent (APPROX) @|, producing rate results that are similar to those 
of HOGWfLD! Recently, Gupta et al. 0 gave an empirical analysis of the effects of a low-precision variant of SGD 
on neural network training. Other variants of stochastic algorithms have been proposed iSB[HEIlmini; only a 
fraction of these algorithms have been analyzed in the asynchronous case. Unfortunately, a new variant of SGD (or a 
related algorithm) may violate the assumptions of existing analysis, and hence there are gaps in our understanding of 
these techniques. 

One approach to filling this gap is to analyze each purpose-built extension from scratch: an entirely new model for 
each type of asynchrony, each type of precision, etc. In a practical sense, this may be unavoidable, but ideally there 
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would be a single technique that could analyze many models. In this vein, we prove a martingale-based result that 
enables us to treat many different extensions as different forms of noise within a unified model. We demonstrate our 
technique with three results: 

1. For the convex case, Hogwild! requires strict sparsity assumptions. Using our techniques, we are able to relax 
these assumptions and still derive convergence rates. Moreover, under HOGWILD !’s stricter assumptions, we 
recover the previous convergence rates. 

2. We derive convergence results for an asynchronous SGD algorithm for a non-convex matrix completion problem. 
We derive the first rates for asynchronous SGD following the recent (synchronous) non-convex SGD work of 
De Sa et al. 0. 

3. We derive convergence rates in the presence of quantization errors such as those introduced by fixed-point 
arithmetic. We validate our results experimentally, and show that Buckwild! can achieve speedups of up to 
2.3x over Hogwild!- based algorithms for logistic regression. 

One can combine these different methods both theoretically and empirically. We begin with our main result, which 
describes our martingale-based approach and our model. 

2 Main Result 

Analyzing asynchronous algorithms is challenging because, unlike in the sequential case where there is a single copy 
of the iterate x, in the asynchronous case each core has a separate copy of x in its own cache. Writes from one core 
may take some time to be propagated to another core’s copy of x, which results in race conditions where stale data is 
used to compute the gradient updates. This difficulty is compounded in the non-convex case, where a series of unlucky 
random events—bad initialization, inauspicious steps, and race conditions—can cause the algorithm to get stuck near 
a saddle point or in a local minimum. 

Broadly, we analyze algorithms that repeatedly update x by running an update step 

x t +i =x t -G t {xt), (2) 

for some i.i.d. update function Gf. For example, for SGD, we would have G(x) = oV/ f (x). The goal of the algorithm 
must be to produce an iterate in some success region S —for example, a ball centered at the optimum x*. For any T, 
after running the algorithm for T timesteps, we say that the algorithm has succeeded if x f G S for some t < T; 
otherwise, we say that the algorithm has failed, and we denote this failure event as Ft. 

Our main result is a technique that allows us to bound the convergence rates of asynchronous SGD and related 
algorithms, even for some non-convex problems. We use martingale methods, which have produced elegant conver¬ 
gence rate results for both convex and some non-convex 0] algorithms. Martingales enable us to model multiple forms 
of error—for example, from stochastic sampling, random initialization, and asynchronous delays—within a single sta¬ 
tistical model. Compared to standard techniques, they also allow us to analyze algorithms that sometimes get stuck, 
which is useful for non-convex problems. Our core contribution is that a martingale-based proof for the convergence 
of a sequential stochastic algorithm can be easily modified to give a convergence rate for an asynchronous version. 

A supermartingale 0 is a stochastic process W t such that E[Wt+i|Wt] < W t . That is, the expected value is 
non-increasing over time. A martingale-based proof of convergence for the sequential version of this algorithm must 
construct a supermartingale W f (xt , x f _ i,..., xo) that is a function of both the time and the current and past iterates; 
this function informally represents how unhappy we are with the current state of the algorithm. Typically, it will have 
the following properties. 

Definition 1. For a stochastic algorithm as described above, a non-negative process Wt : [R rixt —> IR is a rate 
supermartingale with horizon B if the following conditions are true. First, it must be a supermartingale; that is, for 
any sequence Xt,. ■ ■ ,xq and any t < B, 

E[W t+ i{x t - G t (x t ),x t ,... ,x 0 )] < W t {x tl x t - 1 ,... ,a;o). (3) 

Second, for all times T < B and for any sequence xt ,..., xq, if the algorithm has not succeeded by time T (that is, 
Xt £ S for all t < T ), it must hold that 

Wt{xt, xt-i, ■ ■ ■, xq) > T. (4) 

This represents the fact that we are unhappy with running for many iterations without success. 
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Using this, we can easily bound the convergence rate of the sequential version of the algorithm. 

Statement 1 . Assume that we run a sequential stochastic algorithm, for which W is a rate supermartingale. For any 
T < B, the probability that the algorithm has not succeeded by time T is 


P (F t ) < 


E [Wp Qr 0 )] 
T 


Proof. In what follows, we let Wt denote the actual value taken on by the function in a process defined by 0. That 
is, Wt = W t (x t , att-i,..., xo). By applying (0 recursively, for any T, 

E [W T ] < E[W 0 ] = E[W 0 (aro)]. 

By the law of total expectation applied to the failure event Ft, 

E[Wo(* 0 )] > E[Wr] = P(F r )E[W T |F T ] + P(-iF t )E[W t \->F t ]. 

Applying 0, i.e. E[Wt|Et] > 7 , and recalling that W is nonnegative results in 

E[W 0 (x 0 )] >P(F t )T; 

rearranging terms produces the result in Statement|T| □ 

This technique is very general; in subsequent sections we show that rate supermartingales can be constructed for 
SGD on all convex problems and for some algorithms for non-convex problems. 


2.1 Modeling Asynchronicity 

The behavior of an asynchronous SGD algorithm depends both on the problem it is trying to solve and on the hardware 
it is running on. For ease of analysis, we assume that the hardware has the following characteristics. These are basically 
the same assumptions used to prove the original HOGWILD! result & 

• There are multiple threads running iterations of 0, each with their own cache. At any point in time, these 
caches may hold different values for the variable x, and they communicate via some cache coherency protocol. 

• There exists a central store S (typically RAM) at which all writes are serialized. This provides a consistent value 
for the state of the system at any point in real time. 

• If a thread performs a read R of a previously written value X , and then writes another value Y (dependent on 
R), then the write that produced A' will be committed to S before the write that produced Y. 

• Each write from an iteration of 0 is to only a single entry of x and is done using an atomic read-add-write 
instruction. That is, there are no write-after-write races (handling these is possible, but complicates the analysis). 

Notice that, if we let x t denote the value of the vector x in the central store S after t writes have occurred, then since 
the writes are atomic, the value of xt+i is solely dependent on the single thread that produces the write that is serialized 
next in S. If we let Gt denote the update function sample that is used by that thread for that write, and v t denote the 
cached value of x used by that write, then 

x t +i =x t - G t (v t ) (5) 

Our hardware model further constrains the value of vp. all the read elements of v t must have been written to S at some 
time before t. Therefore, for some nonnegative variable T,,t, 


T ~ T 
e i V t = Xt-fi 


( 6 ) 


where e, is the ith standard basis vector. We can think of r i t as the delay in the / 1 h coordinate caused by the parallel 
updates. 

We can conceive of this system as a stochastic process with two sources of randomness: the noisy update function 
samples Gt and the delays fi tt - We assume that the Gt are independent and identically distributed—this is reasonable 
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because they are sampled independently by the updating threads. It would be unreasonable, though, to assume the 
same for the f ? ;. t , since delays may very well be correlated in the system. Instead, we assume that the delays are 
bounded from above by some random variable f. Specifically, if Ft, the filtration, denotes all random events that 
occurred before timestep t, then for any i, t, and k, 

P(n,t>k\Ft)<P(f>k). (7) 

We let r = E[r], and call r the worst-case expected delay. 


2.2 Convergence Rates for Asynchronous SGD 

Now that we are equipped with a stochastic model for the asynchronous SGD algorithm, we show how we can use a 
rate supermartingale to give a convergence rate for asynchronous algorithms. To do this, we need some continuity and 
boundedness assumptions; we collect these into a definition, and then state the theorem. 

Definition 2. An algorithm with rate supermartingale W is (II, R. ^-bounded if the following conditions hold. First, 
W must be Lipschitz continuous in the current iterate with parameter H\ that is, for any t, u, v, and sequence 

Xt , . . . , Xq , 

\\W t (u,x t -i, ...,x 0 )- W t (v,x t - 1 , ■.. ,x 0 )||< H\\u - i>||. (8) 

Second, G must be Lipschitz continuous in expectation with parameter R\ that is, for any it, and v, 

E[||G(«)-G(t;)||]<R||u-t)|| 1 . (9) 

Third, the expected magnitude of the update must be bounded by £. That is, for any x, 

E[||G(aO||]<£. (10) 


Theorem 1 . Assume that we run an asynchronous stochastic algorithm with the above hardware model, for which W 
is a (H , R, £)-bounded rate supermartingale with horizon B. Further assume that II Rfr < 1. For any T < B, the 
probability that the algorithm has not succeeded by time T is 


P (F t ) < 


E[FL(0,xq)] 
(1 - HR^t)T 


Note that this rate depends only on the worst-case expected delay r and not on any other properties of the hardware 
model. Compared to the result of StatementQ] the probability of failure has only increased by a factor of 1 — II Rfr. 
In most practical cases, HR^t <C 1, so this increase in probability is negligible. 

Since the proof of this theorem is simple, but uses non-standard techniques, we outline it here. First, notice that 
the process W t , which was a supermartingale in the sequential case, is not in the asynchronous case because of the 
delayed updates. Our strategy is to use W to produce a new process V* that is a supermartingale in this case. For any 
t and x., if x u ^ S for all u < t, we define 


OO OO 

V t (x t , ...,x 0 ) = W t (x t , ...,x 0 )~ HR^rt + HR^ lkt-fc+i - Xt-k\\ ^ P (t > m). 

k =1 m=k 


Compared with W, there are two additional terms here. The first term is negative, and cancels out some of the 
unhappiness from (|4} that we ascribed to running for many iterations. We can interpret this as us accepting that we 
may need to run for more iterations than in the sequential case. The second term measures the distance between recent 
iterates; we would be unhappy if this becomes large because then the noise from the delayed updates would also be 
large. On the other hand, if x u £ S for some u < t, then we define 


V t (x t , ...,x u ,...,x 0 )= V u (x u ,.. .,x 0 ). 


We call V t a stopped process because its value doesn’t change after success occurs. It is straightforward to show 
that V t is a supermartingale for the asynchronous algorithm. Once we know this, the same logic used in the proof of 
Statement|T|can be used to prove TheoremQ] 
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Theorem |T| gives us a straightforward way of bounding the convergence time of any asynchronous stochastic 
algorithm. First, we find a rate supermartingale for the problem; this is typically no harder than proving sequential 
convergence. Second, we find parameters such that the problem is ( H , R, ^-bounded, typically ; this is easily done 
for well-behaved problems by using differentiation to bound the Lipschitz constants. Third, we apply Theorem|7]to 
get a rate for asynchronous SGD. Using this method, analyzing an asynchronous algorithm is really no more difficult 
than analyzing its sequential analog. 


3 Applications 

Now that we have proved our main result, we turn our attention to applications. We show, for a couple of algo¬ 
rithms, how to construct a rate supermartingale. We demonstrate that doing this allows us to recover known rates for 
Hogwild! algorithms as well as analyze cases where no known rates exist. 


3.1 Convex Case, High Precision Arithmetic 

First, we consider the simple case of using asynchronous SGD to minimize a convex function f(x) using unbiased 
gradient samples V/(.t). That is, we run the update rule 

Xt+l = X t - Q'V ft (x). (11) 

We make the standard assumption that / is strongly convex with parameter c; that is, for all x and y 

(x - y) T (V/0) - V/(y)) > c\\x - yf. (12) 

We also assume continuous differentiability of V/ with 1-norm Lipschitz constant L, 

E[||V/(tr) — V/(y)||] < L\\x — y\\\. (13) 

We require that the second moment of the gradient sample is also bounded for some M > 0 by 

E[||V/>)|| 2 ] < M 2 . (14) 


For some e > 0, we let the success region be 

S = {at|||tc — £*|| 2 < e}- 


Under these conditions, we can construct a rate supermartingale for this algorithm. 

Lemma 1. There exists a Wt where, if the algorithm hasn ’t succeeded by timestep t, 

Wt(x ” ■ ■ ■ ’ X »> = 2ac - ‘° g (‘ 1111 “ ^ £ "‘) + *’ 

such that Wt is a rate submartingale for the above algorithm with horizon B = oo. Furthermore, it is (H,R,f)- 
bounded with parameters: H = 2^/e(2ace — a 2 M 2 ) , R = aL, and £ = aM. 

Using this and Theorem[I]gives us a direct bound on the failure rate of convex Hogwild! SGD. 


Corollary 1. Assume that we run an asynchronous version of the above SGD algorithm, where for some constant 
d € (0,1) we choose step size 

ced 

(\ = - . 

M 2 + 2 LMry/e 

Then for any T, the probability that the algorithm has not succeeded by time T is 


P (F t ) < 


M 2 + 2 LMr^fl 
c 2 edT 



This result is more general than the result in Niu et al. & The main differences are: that we make no assumptions 
about the sparsity structure of the gradient samples; and that our rate depends only on the second moment of G and 
the expected value of f, as opposed to requiring absolute bounds on their magnitude. Under their stricter assumptions, 
the result of Corollary[T]recovers their rate. 
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3.2 Convex Case, Low Precision Arithmetic 

One of the ways BUCKWILD! achieves high performance is by using low-precision fixed-point arithmetic. This 
introduces additional noise to the system in the form of round-off error. We consider this error to be part of the 
BUCKWILD! hardware model. We assume that the round-off error can be modeled by an unbiased rounding function 
operating on the update samples. That is, for some chosen precision factor k, there is a random quantization function 
Q such that, for any at £ R, it holds that E[Q(a;)] = x, and the round-off error is bounded by \Q(x) — a,j< anM. 
Using this function, we can write a low-precision asynchronous update rule for convex SGD as 

x t+ i = x t - Qt (aVft(v t y\, (15) 


where Q t operates only on the single nonzero entry of V ft(vt)- In the same way as we did in the high-precision 
case, we can use these properties to construct a rate supermartingale for the low-precision version of the convex SGD 
algorithm, and then use Theorem[l]to bound the failure rate of convex BUCKWILD! 

Corollary 2. Assume that we run asynchronous low-precision convex SGD, and for some i9 € (0,1), we choose step 
size 

ced 

“ = M 2 ( 1 + k 2 ) + LMt{ 2 + K 2 )ff~e 

then for any T, the probability that the algorithm has not succeeded by time T is 


P (F t ) < 


M 2 ( 1 + k 2 ) + LMt( 2 + n 2 )ffe 
c 2 edT 



Typically, we choose a precision such that k -C 1; in this case, the increased error compared to the result of 
Corollary Q] will be negligible and we will converge in a number of samples that is very similar to the high-precision, 
sequential case. Since each Buckwild! update runs in less time than an equivalent Hogwild! update, this result 
means that an execution of Buckwild! will produce same-quality output in less wall-clock time compared with 
Hogwild! 


3.3 Non-Convex Case, High Precision Arithmetic 


Many machine learning problems are non-convex, but are still solved in practice with SGD. In this section, we show 
that our technique can be adapted to analyze non-convex problems. Unfortunately, there are no general convergence 
results that provide rates for SGD on non-convex problems, so it would be unreasonable to expect a general proof 
of convergence for non-convex Hogwild! Instead, we focus on a particular problem, low-rank least-squares matrix 
completion, 

minimize E[|| A — ata; T ||f.] 
subject to x € IR", 


for which there exists a sequential SGD algorithm with a martingale-based rate that has already been proven. This 
problem arises in general data analysis, subspace tracking, principle component analysis, recommendation systems. 


and other applications 0]. In what follows, we let A = E A . We assume that A is symmetric, and has unit 


eigenvectors ui,U 2 , ■ ■ ■ ,u n with corresponding eigenvalues Ai > A 2 > • • • > A n . We let A, the eigengap , denote 
A = Ai — A 2 . 

De Sa et al. 0] provide a martingale-based rate of convergence for a particular SGD algorithm, Alecton, running 
on this problem. For simplicity, we focus on only the rank-1 version of the problem, and we assume that, at each 
timestep, a single entry of A is used as a sample. Under these conditions, Alecton uses the update rule 


xt +1 = (I + vn 2 e- it e{ t Ae- jt ej t )x t , 


(17) 


where i t and j t are randomly-chosen indices in [1, n\. It initializes x-o uniformly on the sphere of some radius centered 
at the origin. We can equivalently think of this as a stochastic power iteration algorithm. For any e > 0, we define the 
success set S to be 


S = {x\(u^x) 2 > (1 — e) ||x|| 2 }. 


( 18 ) 
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That is, we are only concerned with the direction of x, not its magnitude; this algorithm only recovers the dominant 
eigenvector of A, not its eigenvalue. In order to show convergence for this entrywise sampling scheme, De Sa et al. 
Hi require that the matrix A satisfy a coherence bound m. 

Definition 3. A matrix A £ R" xn is incoherent with parameter p if for every standard basis vector e :j , and for all unit 
eigenvectors Ui of the matrix, (ejut) 2 < /i 2 ro _1 . 

They also require that the step size be set, for some constants 0 < 7 < 1 and 0 < i) < (1 + e ) -1 as 

Aeyd 

V = 2np* \\A\\ 2 f 

For ease of analysis, we add the additional assumptions that our algorithm runs in some bounded space. That is, for 
some constant C, at all times f, 1 < ||xt|| and ||ar* || 1 < C. As in the convex case, by following the martingale-based 
approach of De Sa et al. 0], we are able to generate a rate supermartinagle for this algorithm—to save space, we only 
state its initial value and not the full expression. 

Lemma 2. For the problem above, choose any horizon B such that pyeAB < 1. Then there exists a function Wt such 
that Wt is a rate supermartingale for the above non-convex SGD algorithm with parameters H = 8 nt 7 —1 7 —1 A _ 1 e _ 2 , 
R = pp ||A|| f , and t; = pp ||A|| F C, and 

E [Wo(a:o)] < 2p~ 1 A ~ 1 log(en 7 _ 1 e _1 ) + B^/Ttry. 

Note that the analysis parameter 7 allows us to trade off between B, which determines how long we can run 
the algorithm, and the initial value of the supermartingale E [Wo(a:o)]. We can now produce a corollary about the 
convergence rate by applying Theorem[l]and setting B and T appropriately. 

Corollary 3. Assume that we run HOGWILD! Alecton under these conditions for T timesteps, as defined below. Then 
the probability of failure, P (Ft), will be bounded as below. 


T = 


4 np 4 \\A\\ 2 f 

A 2 e7'i? v / 27r 7 


log 



P (F t ) < 


yffjFyp 2 

p 2 ~\CdTjl' 


The fact that we are able to use our technique to analyze a non-convex algorithm illustrates its generality. Note 
that it is possible to combine our results to analyze asynchronous low-precision non-convex SGD, but the resulting 
formulas are complex, so we do not include them here. 


4 Experiments 


We validate our theoretical results for both asynchronous non-convex matrix completion and BUCKWILD!, a HOG- 
WILD! implementation with lower-precision arithmetic. Like Hogwild!, a BUCKWILD! algorithm has multiple 
threads running an update rule (0 in parallel without locking. Compared with Hogwild!, which uses 32-bit floating 
point numbers to represent input data, Buckwild! uses limited-precision arithmetic by rounding the input data to 
8 -bit or 16-bit integers. This not only decreases the memory usage, but also allows us to take advantage of single- 
instruction-multiple-data (SIMD) instructions for integers on modern CPUs. 

We verified our main claims by running Hogwild! and Buckwild! algorithms on the discussed applications. 
Table Q] shows how the training loss of SGD for logistic regression, a convex problem, varies as the precision is 
changed. We ran SGD with step size a = 0.0001; however, results are similar across a range of step sizes. We 
analyzed all four datasets reported in DimmWitted |25] that favored Hogwild!: Reuters and RCV1, which are 
text classification datasets; Forest, which arises from remote sensing; and Music, which is a music classification 
dataset. We implemented all GLM models reported in DimmWitted, including SVM, Linear Regression, and Logistic 
Regression, and report Logistic Regression because other models have similar performance. The results illustrate that 
there is almost no increase in training loss as the precision is decreased for these problems. We also investigated 
4-bit and 1-bit computation: the former was slower than 8 -bit due to a lack of 4-bit SIMD instructions, and the latter 
discarded too much information to produce good quality results. 
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Table 1: Training loss of SGD as a function of arithmetic precision for logistic regression. 


Dataset 

Rows 

Columns 

Size 

32-bit float 

16-bit int 

8-bit int 

Reuters 

8K 

18K 

1.2GB 

0.5700 

0.5700 

0.5709 

Forest 

58 IK 

54 

0.2GB 

0.6463 

0.6463 

0.6447 

RCV1 

78 IK 

47K 

0.9GB 

0.1888 

0.1888 

0.1879 

Music 

515K 

91 

0.7GB 

0.8785 

0.8785 

0.8781 



threads 


iterations (billions) 


(a) Speedup of BUCKWILD! for dense RCV1 dataset. 


(b) Convergence trajectories for sequential versus HOGWILD! 
Alecton. 


Figure 1: Experiments compare the training loss, performance, and convergence of Hogwild! and Buckwild! 
algorithms with sequential and/or high-precision versions. 


Figure [T(a)| displays the speedup of Buckwild! running on the dense-version of the RCV1 dataset compared to 
both full-precision sequential SGD (left axis) and best-case Hogwild! (right axis). Experiments ran on a machine 
with two Xeon X650 CPUs, each with six hyperthreaded cores, and 24GB of RAM. This plot illustrates that incorpo¬ 
rating low-precision arithmetic into our algorithm allows us to achieve significant speedups over both sequential and 
Hogwild! SGD. (Note that we don’t get full linear speedup because we are bound by the available memory band¬ 
width; beyond this limit, adding additional threads provides no benefits while increasing conflicts and thrashing the 
LI and L2 caches.) This result, combined with the data in Table [Q suggest that by doing low-precision asynchronous 
updates, we can get speedups of up to 2.3 x on these sorts of datasets without a significant increase in error. 

Figure [T(bj| compares the convergence trajectories of Hogwild! and sequential versions of the non-convex Alec¬ 
ton matrix completion algorithm on a synthetic data matrix A £ IR" x " with ten random eigenvalues A, > 0. Each 
plotted series represents a different run of Alecton; the trajectories differ somewhat because of the randomness of the 
algorithm. The plot shows that the sequential and asynchronous versions behave qualitatively similarly, and converge 
to the same noise floor. For this dataset, sequential Alecton took 6.86 seconds to run while 12-thread Hogwild! 
Alecton took 1.39 seconds, a 4.9 x speedup. 


5 Conclusion 

This paper presented a unified theoretical framework for producing results about the convergence rates of asynchronous 
and low-precision random algorithms such as stochastic gradient descent. We showed how a martingale-based rate 
of convergence for a sequential, full-precision algorithm can be easily leveraged to give a rate for an asynchronous, 
low-precision version. We also introduced Buckwild!, a strategy for SGD that is able to take advantage of modern 
hardware resources for both task and data parallelism, and showed that it achieves near linear parallel speedup over 
sequential algorithms. 
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A Proof of Theorem U 

Proof of The orem\J\ This proof is a more detailed version of the argument outlined in Section [2721 First, we restate 
the definition of the process V t from the body of the paper. As long as the algorithm hasn’t succeeded yet. 


771=1 


V t {x t , ...,x 0 ) = Wt{x t , ■ ■ • ,x 0 ) - HR^rt + HR^_ \\x t -k+i - x t -k\\ Y P (f > m). 

k =1 m=k 

At the next timestep, we will have Xt+i = Xt + G(vt), and so 

OO 

V t+ i{x t + G(v t ),x t ,... ,x 0 ) = W t +i(x t + G(v t ), x t , ...,x 0 )~ HR£r(t + 1) + HR G(t5 t )| Y, P (f > to) 

OO OO 

+ HR E ||*£t— fc+2 %t— fc+l|| E P (f > m). 

k—2 m=k 

Re-indexing the second sum and applying the definition of r produces 

V t+ i {x t + G(v t ), x t , ...,x 0 )= W t +i{x t + G(v t ),x t , ...,x 0 )~ HR£r(t + 1) + HRt G(v t ) 

OO OO 

+ HR E fc+l %t— k\\ E P (f > to) . 

k= 1 m=k-\-l 

Applying the Lipschitz continuity assumption © for W results in 

V t +i (x t + G(v t ), x t , ..., xo )< W t+ i(x t + G(x t ), x t , ...,x 0 ) + H G(v t ) - G(x t ) - HR£r(t + 1) 

+ HRt 


G{v t ) +HR ||x t _ fc+ i - x t -k || Y P (t > m). 

m=k-\-l 


k=l 


Taking the expected value of both sides produces 


E 


V t+ i(x t + G(v t ),x t ,.. .,x 0 ) 


< E 


W t+ i(x t + G(x t ), x t ,...,x 0 ) 


H E 


G(v t )-G(x t ) — HR^r(t + 1) 


+ HRtE 


G(v t ) + HR E \\ x t-k+i — x t ~k\\ Y P^>m). 

m=k -\-1 


k=l 


Applying the rate supermartingale property © of W, 


E 


Vt(xt + G(v t ),xt,...,xo) < W t (x t ,... ,x 0 ) + HE G(v t ) — G(x t ) -HR£r(t + 1) 

+ HRtE 


G(v t ) J + HR Ell*t-fc+i-*t-fc|| E P ^>m). 

m=k-\-l 


k=1 


Applying the Lipschitz continuity assumption © for G, 


E 


V t {x t +G{v t ),x t ,...,x 0 ) < W t {x t ,...,x 0 ) +HRE[\\v t - x t \IJ - HRfr(t + 1) 


HRt E 


G(v t ) 1 + HR Y\\xt-k+i-x t -k\\ E P (T>m). 


k= 1 


m=k-\-l 


li 






















































Finally, applying the update distance bound ( fTTTb . 


E 


Vt{x t + G(v t ),x t , ...,x 0 ) < W t (x t ,... ,x 0 ) + HKE [||6 t - XtUJ - HR£r(t + 1) 


HR£t + HRJ2 


\ x t-k -\-1 %t—k\ 


E P(f>m) 


k =1 


m=k -\-1 


= W t (x t ,... ,x 0 ) - HR£rt + HR^ \\xt-k+i ~ x t -k\\ Yf P (t > m) 

k— 1 m—k 

oo 

+ HKE [\\v t - EtHJ - HRJ2 \\ x t-k+i ~ x t-k\\ P (r > k) 

k=l 

oo 

= V t (x u . .. ,xq) + HKE [|| v t ~ II x t-k+i ~ x t -k\\ P(r>k). 


k =l 


Now, by the definition of the v t . 


\vt-xt\h 


= E 

i=1 


T T~ I 

e i x t ~ e i v t \ 


, T I 

| e i x t e i X t—Ti, t | 


= E 

i= 1 
n Ti,t 

< E E \ e i x t-k+i - ejx t - k 


2— 1 /C = l 

Furthermore, using the bound on f^t from (0! gives us 


n oo 

ef xt-k+i - ef x t -k | P {n, t > k ) 

2—1 /C—1 

n oo 

^EEI ef a; t _ fe+ i - ef x t _ fc | P (f > k) 

2—1 k—1 
oo 

= ^2 \\ x t-k+i - x t-k 111 P (f > k) 
k= 1 
oo 

= ^ ||aj t — fc _|_i - ar t _fc|| P (f > k), 

fc=i 


where the 1-norm is equal to the 2-norm here because each step only updates a single entry of x. Substituting this 
result in to the above equation allows us to conclude that, if the algorithm hasn’t succeeded by time t. 


E 


V t (x t + G(v t ),x t , ...,x 0 ) 


< Vt(x t , ■ ■ -,x 0 ). 


(19) 


On the other hand, if it has succeeded, this statement will be vacuously true, since Vt does not change after success 
occurs. Therefore, (fl9] > will hold for all times. 

In what follows, as in the proof of Statement 0 we let V t denote the actual value taken on by the function during 
execution of the algorithm. That is, Vt = V t (xt,Xt- 1 ,..., Xq). By applying (fl~9] > recursively, for any T < B, we can 
show that 

E [V T ] < E [V 0 ] ■ 

Since we assumed as part of our hardware model that Xt = Xq for t < 0, 


E [Vo] = E [Wb(a:o)]. 
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Therefore, by the law of total expectation 


E [W„0r o )] 

> E [V T \ 

= E [V t \F t ] P (F t ) + E [V T \->F T \ P (~>F T ) 

> E[E T | Ft)P(Ft) 


= E 


W T (x T , ...,x 0 )- HR£tT + HR E \\XT-k+l - %T—k || ^ P (f > m) 


k=l 


m=k 


Ft 


P (F t ) 


> (E [W T (x T , • • •, x 0 )\Ft\ ~ HR^tT) P (F t ) . 


Since Wt is a rate supermartingale, we can apply (Q} to get 

E [Wo(i 0 )] > (T - HR£tT) P (F t ) , 


and solving for P (Ft) produces 


as desired. 


P (Ft) < 


E[W 0 (a;o)] 

(1 - HR^t)T ’ 


□ 


B Proofs for Convex Case 

First, we state the rate supermartingale lemma for the low-precision convex SGD algorithm. 
Lemma 3. There exists a Wt with 


Wo(zo) < 


e 

2 ace — a 2 M 2 ( 1 + n 2 ) 


log 


||tr 0 -x*|| 2 ^ 


such that Wt is a rate submartingale for the above convex SGD algorithm with horizon B = oo. Furthermore, it is 
(H, R, £)-bounded with parameters: R = aL, £ 2 = a 2 ( 1 + n 2 )M 2 , and 


H = 


2 yfe 

2ace — a 2 M 2 ( 1 + n 2 )' 


We note that, including this Lemma, the results in Section l3Tl are the same as the results in Section [L2l except that 
the quantization factor is set as k = 0. It follows that it is sufficient to prove only the Lemma and Corollary in 13.21 
this is what we will do here. 

In order to prove the results in this section, we will need some definitions and lemmas, which we state now. 


Definition 4 (Piecewise Logarithm). For the purposes of this document, we define the piecewise logarithm function 
to be 

... / log(ex) : x > 1 

l0gW = \ * :*<1 

Lemma 4. The piecewise logarithm function is differentiable and concave. Also, if x > 1, then for any A, 


log(a;(l + A)) < log(x) + A. 

Proof The first part of the lemma follows from the fact that log (a;) is a piecewise function, where the pieces are both 
increasing and concave, and the fact that the function is differentiable at x = 1. The second part of the lemma follows 
from the fact that a first-order approximation always overestimates a concave function. □ 

Armed with this definition, we prove Lemma[3] 
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Proof of Lemma\3\ First, we note that, at any timestep t , if we evaluate the distance to the optimum at the next timestep 
using CD, then 

2 


Xt + G t (x t ) - x* 

= II x t - x*|| 2 - 2 (x t - x*) T Q t (aV/t(x t )) + 

Qt (aVMxtj) 


= ||x t - x*|| 2 - 2(x t - x*) T Qt ( aVf t (x t )) 



aV/ t ( x t ) + Qt (c*V ft(xf - aV f t (x t ) 


Taking the expected value and applying (fl4l >. and the bounds on the properties of Qt , produces 


E 


x t + G t (x t ) - x* 


< || x t — x*|| — 2 a(xt — X *) 1 V/(xt) + a 2 M 2 + S 2 . 


Since we assigned 6 < anM, 


E 


x t + G t (x t ) - x* 


< II Xt - x*|| 2 - 2 a(x t - x*) T S7f(x t ) + a 2 M 2 ( 1 + n 2 ) 

= || x t - x*|| 2 - 2a(x t - x*) T (V/(® t ) - V/(®*)) + a 2 M 2 { 1 + n 2 ). 

Applying the strong convexity assumption (fl2] >. 

< ||Xt — tr*|| 2 — 2ac||x t — x*|| 2 + a 2 M 2 ( 1 + n 2 ) 

= (1 — 2 ac) ||xt — x*|| 2 + a 2 M 2 ( 1 + n 2 ). 

Now, if we haven’t succeeded yet, then ||x t — x*|| 2 > e. Under these conditions, 

2 


E 


x t + G t (x t ) - x* 


E 


x t + G t {x t ) - x* 


< ||xt — x*|| (l — 2ac + a 2 M 2 (l + n 2 )e 1 ) . 


Multiplying both sides of the equation by e 1 and taking the piecewise logarithm, by Jensen’s inequality 


E 


log e 


-1 


x t + G t (x t ) - x* 


< log E 


,-i 


x t + G t (x t ) - x* 


< 


log (e 1 ||x t — x*|| 2 (l — 2ac + a 2 M 2 (\ + n 2 )t x )) • 


Since e 1 ||x t — x*|| 2 > 1, we can apply LemmaQ] which gives i 


E 


log e 


t 


x t + G t (x t ) - x* 


< log (e 1 ||x t — x*|| 2 ^ — 2ac + a 2 M 2 (l + n 2 )e 1 . 


Now, we define the rate supermartingale Wt such that, if we haven’t succeeded up to time t, then 


W t {x t , ...,x 0 ) = 


■ log (e Ml 


x t - x 


2 ace — a 2 M 2 ( 1 + k 2 ) 
otherwise, if u is a time such that x u £ S, then for all t > u, 

W t (x t , ■ •• ,x 0 ) = W u {x u , ■ ■ • ,x 0 ). 

The first rate supermartingale property © is true because if success hasn’t occurred, 

e 


+ t'i 


E 


Wt+i(x t + Gtixt ),..., x’o) 


= E 


< 


2ace — a 2 M 2 {l + k 2 ) 
e 


E 


j log 


x t + G t (x t ) - x* 

log 1 

v. 

x t + G t (x t ) - x* 


+ {t + 1) 


(log ( e 1 \\xt~x* II 2 ) 


2ace — a 2 M 2 { 1 + k 2 ) 

+ (t + 1) 

= 2a,„ - aHP( 1 + «*) l»g l|x, - ) - 1 + (t + i) 

= W t (x t ,... ,x 0 ); 


— 2ac+ a 2 M 2 (l + n 2 )e~ 1 ^j 
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it is vacuously true if success has occurred because the value of Wt does not change after x u £ S for u < t. The 
second rate supermartingale property IQ} holds because, if success hasn’t occurred by time T, 


W T (x T , ■ ■ ■ ,x 0 ) = 


2 ace — a 2 M 2 ( 1 + k 2 ) 


log (e ^lIxT-aJ*!! 2 ) 


+ T > T; 


this follows from the non-negativity of the log function for non-negative arguments. 

We have now shown that Wt is a rate supermartingale for this algorithm. Next, we verify that the bound on Wq 
given in the lemma statement holds. At time 0, by the definition of the log function, since we assume that success has 
not occurred yet. 


Wo(x 0 ) = 


2 ace — a 2 M 2 ( 1 + k 2 ) 
e 

2 ace — a 2 M 2 ( 1 + k 2 ) 


log(e 1 ||x 0 -x*|| 2 ) 


log 


e ||x 0 - x*||" 


this is the bound given in the lemma statement. 

Next, we show that this rate supermartingale is (H,R,£)- bounded, for the values of H , i?, and £ given in the 
lemma statement. First, for any x, t, and sequence x t - \ ..... .x'o. 


V x W t (x,x t - 1 ,... ,x 0 ) = Vj 


2ace — a 2 M 2 ( 1 + n 2 ) 

a 

- 26 - 


log (e 1 ||x-x*|| 2 )^ 

1 (x x*)log ^e _1 1|x — x*|| 2 ^ . 


2ace — a 2 M 2 { 1 + k 2 ) 

Now, by the definition of log, we can conclude that log 7 (it) = min (l, it -1 ). Therefore, 

2 


V x W t (x,x t - 1 ,..,,x 0 ) = 

and taking the norm of both sides, 

V x Wt(x,xt-i, ■ • ■ ,x 0 ) = 


2ace — a 2 M 2 ( 1 + k 2 ) 


2ace — a 2 M 2 (l + k 2 ) 

Clearly, this expression is maximized when ||x — x* || 2 = e. Therefore, 

2\fk 


(x — x*) min ^l,e ||x — x*|| 2 ^j , 


in (||x — x*|| ,e ||x — x* || . 


V x W t (x,x t _i,... ,x 0 ) < 


2ace — a 2 M 2 (l + k 2 ) 


The Lipschitz continuity expression with H in the lemma statement now follows from the mean value theorem. 

Next, we bound the Lipschitz continuity expression for R. We have that, for any x and y, if the single non-zero 
entry of V/ is at index i, then 


E 


G(x) - G(y) 


= E 
= E 


Q(aV/(x)) - Q(aV/(y)) 
Q(aefV/(x)) - Q(aeJS7f(y)) 


Without loss of generality, we assume that Q is non-decreasing, and that ef V/(x) > ej\7 f{y). Thus, by the unbiased 
quality of Q , 


E 


G(x) - G{y) 


= E 
= E 

= aE 


Q(efaVf{x)) - Q(efaVf(y)) 
efaV/(x) - efaVf(y) 

V/(x) - V/(y) 


Finally, applying (IT3l) . 


E 


G(x) - G(y) 


< aL. 
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Finally, we bound the update expression with £. We have, 

2 

E 


G(x) 

= E 
< E 


= E 


Q(aV/(x)) 


Q(aVf{x)) 


V/>)|| 2 + 2a(V/(x)) T (Q(aV/>)) - aV/>)) + ||Q(aV/(x)) - aV/>) 


Applying the bounds on the rounding error, 


E 


G(a 


< E 


V/>)|| 2 + 2a(V/(x)) T (Q(aV/>)) - aV/(x)) + 5 2 


Taking the expected value and applying (fl4l > and the unbiased quality of Q, 

1 2 


E 


Applying the assignment <5 = anM results in 


E 


G(x) 


< a 2 M 2 + S 2 . 


G{ x) <a 2 M 2 { 1 + k 2 ) 


which is the desired expression. 

So, we have proved all the statements in the lemma. 

Proof of Corollary^ Applying TheoremQ]directly to the result of Lemma[T|produces 

E[W 0 (x 0 )] 


□ 


P (Ft) < 


(1 - HR£t)T 
e 

2 ace — a 2 M 2 ( 1 + ft 2 ) 


log 


e ||a? 0 - x 


*||2 


2i/e 


2ace — a 2 M 2 (l + k 2 ) 


< 


(2 ace — a 2 (M 2 ( 1 + ft 2 ) — 2 LMt\J\ + n 2 y / e)) T 
e 


log 


(aL)(aM \/l + k 2 )t ) T 
e ||x 0 - x*|| 2 \ 


(2ace — a 2 (M 2 (l + ft 2 ) — LMt( 2 + ft 2 )i/e)) T 

Substituting the chosen value of a, 


log 


e ||x 0 - x*||" 


P(F T )<-\2ce 


cet? 


T ^ V^ 2 (l + k 2 ) + LMt (2 + ft 2 )^ 

- (M 2 (l + ft 2 ) - LMt(2 + k 2 )v^) ^ 


cet? 


M 2 (l + ft 2 ) + LMt( 2 + n 2 )yfe 


* m2 


log 


e ||x 0 - x* || 


( 


2c 2 e 2r d 


c 2 e 2r d 2 


M 2 (l+K 2 )+LMr(2+K 2 )Ve M 2 (1 +k 2 )+LMt(2+. 

2 




T 


■log 


e ||x 0 - x*||" 


< 


c 2 e 2 i9 


T 


log 


e 11x 0 - x*||" 


M 2 (1+k 2 )+LMt(2+k 2 ) v /£ 

M 2 ( 1 + K 2 ) + LMt(2 + n 2 )y/e f e ||x 0 - x*|| 2 
^&T ° S e 


as desired. 


□ 
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C Proofs for Non-Convex Case 


In order to accomplish this proof, we make use of some definitions and lemmas that appear in De Sa et al. 10. We 
state them here before proceeding to the proof. 

First, we define a function 


T ( x ) = 


(ujx) 2 


(1 — 7 n - 1 )(W[ x) 2 + 7 n -1 


Clearly, 0 < r(.x') < 1. Using this function, De Sa et al. 0 prove the following lemma. While their version of the 
lemma applies to higher-rank problems and multiple distributions, we state here a version that is specialized for the 
rank-1, entrywise sampling case we study in this paper. (This is a combination of Lemma 2 and Lemma 12 from De Sa 
et al. 0 .) 


Lemma 5 (r-bound). If we run the Alecton update rule using entrywise sampling under the conditions in Section ITT 
including the incoherence and step size assignment, then for any x ^ S, 

E t(x + pAx) > t(x) (1 + rjA(l — . 

We also use another lemma from De Sa et al. 0. This is a combination of their Lemmas 1 and 7. 

Lemma 6 (Expected value of r(a;o)). If we initialize Xo with a uniform random angle (as done in Alecton), then 

E [1 — t(xq)\ < 

Now, we prove Lemma[2] 

Proof of Lemma\2\ First, if x ^ S, then (ufx) 2 < (1 — e) ||tc|| 2 . Therefore, 

(ufx) 2 


-(x) = 


< 


(1 — 7 n~ 1 )(uj'x) 2 + 7 n -1 ||a:|| 
1 - e 

(1 — 7 n -1 )(l — e) + 7 n~ l 
1 - e 


and so 


1 — e + 7 n 1 e 
1 — r( x ) > 


7 n 1 e 


1 — e + 7 ?t, 1 e. 


> 771 1 e. 


From the result of Lemma |5] for any x f S, 


E 


t(x + pAx) > t(x) (1 + ? 7 A (1 — t(x))) . 


Therefore, 


E 


1 — t ( x + r/Ax) < (1 — t(x)) (1 — pAr(x)) 


Therefore, by lensen’s inequality and Lemma|4] since 7 1 ne( 1 — t(x)) > 1 , 


E 


log (7 1 ne 1 ^1 — t(x + rjAx)^ > log ^E 7 1 ne 1 ^1 — t(x + tjAx)^ 'j 

> log ( 7 - 1 ne -1 (l — t(x)) (1 — t]At(x))) 

> log (7 _1 7ie _1 (l — t(x))) — 7|At(i). 


Now, we define our rate supermartingale. First, define 


Z = { x 


r (z) > 
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and let B > 0 be any constant. Let W t be defined such that, if x u ^ S U Z for all u < t, then 
W t (x t ,...,x 0 ) = ^-log (7 ~ 1 ne~ 1 (l- T(x t ))) + 2B(1 — r(x t )) + t. 
On the other hand, if x u G S' U Z for some u , then for all t > u, we define 


W t (x t , ...,x 0 ) = W u (x u ,... ,Zo). 


That is, once xt enters S J Z, the process W stops changing. 

We verify that Wt is a rate supermartingale. First, ([3]) is true because, in the case that the process has stopped it is 
true vacuously, and in the case that it hasn’t stopped (i.e. Xi ^ S U Z for all u < t). 


E 


W t+ i ( x t + r/AtXt ,x t ,...,xo 


= E 


log (7 1 ne 1 (l-T(x t +V^tXt))^+‘2B(l-T{x t +r]AtXt)) 


+ t + 1 


= 4 ^ 

?/A 


log ^7 lne 1 (1~ r(x t + rjA t Xt))^ + 2BE 1 - r(x t + rjA t x t ) 


+ t 


+ 1< — (log (7 1 ne 1 {l-T(xt)))-r)AT(xt))+2B(l-T(xt)) + t+l = Wt(x t ,...,Xo)-2T(xt) + l. 


Since x t 0 Z, it follows that 2r(x t ) > 1. Therefore, 


E 


W t +i{x t +r]A t xt,x t , ...,x 0 


< W t (x t , ■. .,X 0 ). 


And so (|3} holds in all cases. 

The second rate supermartingale property (|4} holds because, if success hasn’t occurred by time T < B, then there 
are two possibilities: either the process hasn’t stopped yet, or it stopped at a timestep where Xt G Z. In the former 
case, by the non-negativity of the log function, 

W T {x T , ■ ■ ■, xo) = ^ log ( 7 -1 ne -1 (l - t(x t ))) + 2B(1 - t(x t )) + T >T. 

In the latter case, 

W T {x T , ...,x 0 ) = ^log - t(x t ))) + 25(1 - t(x t )) + T 

> B. 

Therefore © holds. 

We have now shown that Wt is a rate supermartingale for Alecton. Next, we show that our bound on the initial 
value of the supermartingale holds. At time 0, 

W 0 (x 0 ) = log ( 7 -1 ne -1 (l - t(x 0 ))) + 25(1 - t(x 0 )) 

< ^ log ( 7 -1 ne -1 ) + 25(1 - t(x 0 )) 

= ^ log + 2-8(1 - t(x 0 )). 

Therefore, applying LemmaO 

E [W 0 (a;o)] < ^ log + 25E [1 - r(x 0 )] 


This is the value given in the lemma. 
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Now, we show that W t is (//, 11. ^-bounded. First, we give the II bound. To do so, we first differentiate t(x). 


Vr(cc) 


2u\ufx ^(1 — 'yn 1 )(m^’;c ) 2 + jn 1 ||x|| 2 ^ — 2 (w^a ;) 2 ((1 — 7 ?r 1 )uiufx + jn 1 x) 

((1 - 7 rc _ 1 )(wfa ;) 2 + 7 n -1 ||x|| 2 ^ 

2uiujx'yn~ 1 ||x || 2 — 2 (M^cc) 2 7 n _ 1 a; 

((1 - 7 rc _ 1 )(wf;r ) 2 + 7??. -1 ||x|| 2 ^ 

_! uiujx ||x|| 2 - x(u[x) 2 
Z'yn - 2 • 

^(1 — 7 n^ 1 )(rtfa ;) 2 + 7 n ^ 1 ||at||"^ 


Therefore, 

||vr(,)|| 2 = 472n - 2 ^WN^M^NL 
^(1 - 7n -1 )(wfx) 2 + 7n _1 ||x|| 2 ^ 

_ 2 -2 |N| 4 -«x) 2 ||a:|| 2 

< 47 n - 

((1 - 7n _1 )(M^x) 2 + 7 nr 1 ||x|| 2 ^ 

< 4 7 „->--- 

((1 - 7 Jj _ 1 )(wf 7 ) 2 + 7 n _1 ||7|| 2 t 

< _ 4(1 -t(s)) _ 

((1 - 7n _1 )(ufat) 2 + 7n -1 ||x|| 2 ^ 

4 ?r(l — r(at)) 

7 INI 2 


(uJxYWxt-julxYWxW 2 
^(1 - 7 n _ 1 )(wf;r ) 2 + 7 n -1 ||at|| 2 ^ 

iixr-c^iNi 2 

((1 - 7 n _ 1 )(M^x ) 2 + 7 n _1 ||at|| 2 ^ 

ll^ll 2 (1 — r(£c)) 


Applying the assumption that INI > 1 . 


||Vr(a;)|| < 


1 4?r(l — t(x)) 


Now, differentiating Wt with respect to r produces 


dW 2n ' , , 

17 = ’^ bg ( 7 -«- (■ l~r))~ 2 B. 


So, it follows that 

||V x Wt(a;, xt~i, ■. ■, cco)|| < 

< 


dW 


dr 

( 2 n 


IIV -r(at) || 


\77eA l0g (^” lne_1 ( 1 “ T )) + 23 


1 4n(l — t(x)) 


Applying our assumption that r]"/eAB < 1, it is clear that this function will be maximized when 7 1 ne 1 (1 —r) = 1. 
Therefore, 

/ 2n \ 

||V x Wt(at, xt-i, ..., tc 0 ) || < + 2B ) 2 v / e 

8 n 
777 A 


which is our given value for H. 
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Next, we give the R bound. For Alecton, we have 


Therefore, 


E 


G{x) = rjAx = rjnr'eiel Ae^e^ x. 


G(x) - G{y) = ?7?r 2 E [||aef AejeJ (x - y)\\] 
= ?y?r 2 E [| ef AejeJ(x — y)|] 

n n 

= vJ2'52\ef A e j \ \ej(x — y)\ 

*= 1 3= 1 




3 =1 


vi=l 


< l e J(^ — S/)| 


3 =1 




= v 1 _ y) I ^ ( e J ^ 2e J) ’ 


j=i 


= »?X) l e J( a;_ y)l V" ( X) A i( u fe e j ) 2 

f=l 


v.k=l 


Applying the incoherence bound, 
E 


G(z)-G(y) <77^|eJ(x-y)|\/^ ^A 2 /r 2 n 1 


j=i 


V/c—1 


= 1 II^H ' 


. 7=1 


= r7El^-y)|H|A||, 

f=i 

= 77/x||A|| F ||a:-y|| 1 . 


This agrees with our assignment of f? = 77 /r ||^n F . 

Finally, we give our £ bound on the magnitude of the updates. By the same argument as above, we will have 


E 


G(x) = ? 7 n 2 E [||e,;ef AejeJa:||] 
= riH\\A\\ F \\x\\ 1 . 


Applying the assumption that < C, produces the bound given in the lemma, £ = 7777 ||^n F 


C. 


This completes the proof of the lemma. 

Next, we prove the corollary that gives a bound on the failure probability of asynchronous Alecton. 
Proof of Corollary\3\ By TheoremQ] we know that for the constants defined in Lemma[2] 

P(F T )< mAnL . 
k ’ - (1 - HR(r)T 

If we choose B = T for the horizon in Lemma[2] and substitute in the given constants, 

jA log (^) + T ^n) (l - (^7;) (W ll-lll f) (W Mil rC) r) T-> 


□ 


p (Ft) < 
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Now, for the given value of r/, we will have 


8pnp 2 \\A\\ f Ct Aeyd 8npt 2 \\A\\ 2 f Ct 


2||A|| 2 f jA sfe 

4 Ctiryfe 


Also, for the given values of // and T, we will have 

2 


/ en\ 2 np, 4 ||A|| 2 F A 2 e r yd^/2ny 2 

i]AT °^\7 e / Aejd 4np 4 m|| 2 F A 

= \/ 27r 7 - 


Substituting these results in produces 


p(FT)< v ^(i-^AA') 

a / 87 TJfl 2 


-1 


/i 2 - 4Ci?rVe’ 


which is the desired result. 


□ 


D Simplified Convex Result 

In this section, we provide a simplified proof for a result similar to our main result that only works in the convex case. 
This proof does not use any martingale results, and can therefore be considered more elementary than the proofs given 
above; however, it does not generalize to the non-convex case. 

Theorem 2. Under the conditions given in Section ED for any e > 0, if for some d £ (0,1) we choose constant step 
size 


cde 


then there exists a timestep 


such that 


T < 


2LMtJI + M 2 ’ 
2LMr^/e + M 2 1 ( ||cco — ar 


c 2 de 


■log 


E 


\Xt — x 


< e. 


Proof Our goal is to bound the square-distance to the optimum by showing that it generally decreases at each timestep. 
We can show algebraically that 

im+i - x*\\ 2 = im - a;*|| 2 - 2 a{x t - x*) T Vf t (x t ) + 2a(x t - x*) T (\7 f t (x t ) - V/ t (ui)) + « 2 V/ t (u t ) 

We can think of these terms as representing respectively: the current square-distance, the first-order change, the noise 
due to delayed updates, and the noise due to random sampling. Taking the expected value given v t and applying 
Cauchy-Schwarz, ( fl2l . (fl3l> . and (fl4l > produces 


E 


Ft+i - x 


F t ,Vt 


< \\x t - x*\Y - 2ac\\x t - x*|| 2 + 2aL \\x t - tc*|| ||a: t - u t |li + a 2 M 2 

n 

= (1 - 2 ac) \\x t - x*|| 2 + a 2 M 2 + 2 aL \\x t - £*|| ^ \e[x t - ejx t -f it , 


2—1 
n th 


< (1 - 2 ac) || x t - x*|| 2 + a 2 M 2 + 2aL \\x t - *iEE | %t—k -\-1 %t—k \ • 


2—1 k= 1 
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We can now take the full expected value given the filtration, which produces 


E Hxt+i — x*|| 2 Tt < (l — 2ac) \\xt~ x*\\ 2 + a 2 M 2 + 2aL \\xt — x*\\ EE P(fi,k > k) | ej xt-k+i — e[ x t ~k \ ■ 

i= 1 k=l 

Applying 0 results in 

n oo 

E \\x t +i - x*\\ 2 T t < (1 - 2 ac) \\x t - x*\\ 2 + a 2 M 2 + 2aL \\x t - **11 EE^^ - \ e f x t-k+i - ef x t -k\ 

2—1 k =1 
oo 

= (1 - 2ac) II x t - x*\\ 2 + a 2 M 2 + 2 aL \\x t - z*|| ^ P (f > k) ||x t _ fc +i - x t -k\\ 1 , 

and since only at most one entry of x changes at each iteration. 


fc=l 


E 


\\x t +i - x*|| 2 Ft < (1 - 2ac) ||x t - x*|| 2 + a 2 M 2 + 2aL^ P(f > fc) ||a; t - z*|| ||x t _ fc+ i - . 


fe=i 


Finally, taking the full expected value, and applying Cauchy-Schwarz again. 


E 


\xt+i — x*\\ 2 < (1 - 2ac)E p|x t - x*\Y 


+ c?M 2 


Noticing that, from (fl4l >. 


E 


fe=l 


\Xt-k-\-l Xt—k | 


/— 




V E 

II Xt - z*|| 2 

E 

k-\-l Xt— fc|| 


= E 


aG(v t -k ) 


< a 2 M 1 


if we let Jt = E 


n*t-**ir 


, we get 


,/ t+ i < (1 - 2ac)Jt + a 2 M 2 + 2a 2 LM s £ j P(f > k) \J7 t 

k= 1 

= (1 - 2ctc) J t + a 2 M 2 + 2a 2 LMTy/7 t . 

For any e > 0, as long as J t > e, 

log Jt- i-i < log J t + log ^1 — 2ac + a 2 M 2 e ~ 1 + 2 a 2 LMre~^ 
< log J t — 2ac + a 2 M 2 e ~ 1 + 2 a 2 LMre~i. 

If we substitute the value of a chosen in the theorem statement, then 

c 2 t?e 


log J t +1 < log J t - 


2LMt^I + M 2 ' 


Therefore, for any T, if Jt > e for all t <T, 


„ 2LMt+ M 2 f J 0 
1 < -—-log ' 


c 2 tie 


Jt 


which proves the theorem. 


□ 
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